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Abstract 

By numerically calculating the time-evolved Wigner functions, we investigate the dynamics of a few- 
photon supeiposed (e.g., up to two ones) state in a dissipating cavity. It is shown that, the negativity 
of the Wigner function of the photonic state unquestionably vanishes with the cavity's dissipation. As a 

cr 

consequence, the nonclassical effects related to the negativity of the Wigner function should be weakened 
gradually. However, it is found that the value of the second-order correlation function g^(0) (which 
^ ' serves usually as the standard criterion of a typical nonclassical effect, i.e., g^(0) < 1 implies that the 

CO ■ photon is anti-bunching) is a dynamical invariant during the dissipative process of the cavity. This feature 

is also proven analytically and suggests that g^ (0) might not be a good physical parameter to describe the 
photonic decays. Alternatively, we find that the anti-normal-order correlation function g^ 2A ^ (0) changes 
with the cavity's dissipation and thus is more suitable to describe the dissipative-dependent cavity. Finally, 
^ ■ we propose an experimental approach to test the above arguments with a practically-existing cavity QED 

system. 
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I. INTRODUCTION 



It is well-known that the Wigner function, introduced 70 years ago by Wigner to describe the 
quasi-probability distribution of a quantum particle in its phase space, is a very popular tool to 
study the statistical properties of various quantum states [1]. Basically, once the Wigner function 
has been determined, all the knowable information on the quantum state (such as its nonclassical 
statistical properties) can be extracted [2-5]. Typically, differing from the standard probability dis- 
tribution, such a quasi-probability distribution can be assigned by a negative value [6]. Therefore, a 
quantum state with the negative Wigner function should be nonclassical and thus certain nonclassi- 
cal effects (such as the photon anti-buchings) [7-11] demonstrate. This indicates that, determining 
the Wigner function of a selected quantum state plays an important role both fundamentally and 
practically in quantum-state engineerings. 

Usually, any selected quantum system is always surrounded by the classical environments. 
Thus, dissipation of the artificially-prepared quantum sate is one of the central topics in quantum 
coherence science. Roughly, due to the existence of various dissipations and fluctuations from the 
environments, any excited quantum state will decay to the ground state and the relevant system 
finally becomes classical. Under the standard logic, people pay the most attention to calculate 
either decoherence or decay time of a superposition quantum state, rather than cares on the process 
of the decoherence or decay [12-15]. Alternatively, in the present work we investigate exactly the 
dissipative dynamics for a prepared quantum state by calculating its dissipative-dependent Wigner 
function. Our discussions are based on the typical few-photon quantum state in a cavity, but can 
be directly generalized to other quantum systems such as qubits and qutrits. 

The paper is organized as: in Sec. 2, we describe how the Wigner function for a few-photon 
superposed state changes with the cavity's dissipation. Our numerical results show naturally that 
the negativity of the Wigner function weakens gradually with the dissipation and the final state of 
the cavity should be "classical" with positive Wigner function. With the calculated Wigner func- 
tion we investigate how the nonclassical properties, such as the anti-bunching effect of photons, 
changes with the cavity dissipation. It is surprised that the value of the second-order correla- 
tion function g^(0) (which serves usually as the standard criterion of a nonclassical effect, i.e., 
</ 2 )(0) < 1 implies that the photon is anti-bunching) is an invariant during the dissipative pro- 
cess of the cavity. We prove such an argument analytically by directly solving the relevant master 
equation and suggests that g^ 2 \0) is not a good parameter to describe dissipative-dependent non- 
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classicality of the photonic decays. Alternatively, we find that the anti-normal-order correlation 
function g^ 2A \0) changes with the cavity's dissipation and thus could be more suitable to de- 
scribe the dissipative-dependent cavity. With an experimentally-demonstrated cavity QED system 
we propose an approach to test our results, including how to prepare the investigated few-photon 
superposed state of the cavity and measure its Wigner function. Finally, our conclusions and 
discussions are given in Sec. 4. 

II. DISSIPATIVE DYNAMICS OF WIGNER FUNCTIONS FOR FEW-PHOTONS SUPERPOSI- 
TION STATES 

Generally, the quasi-probability distribution W(a, a*) can be defined by the Fourier transform 
of the symmetrical-ordered characteristic function C(A, A*) [16], i.e., 

W(a,a*) = J d 2 \ C(\,\*)e aX *- a * x , (1) 

with A and a being the complex parameters in phase space. The expression of the symmetrical- 
ordered characteristic function is defined as 

C(A,A*) =Tr[pe Aat ~ aA *], (2) 

where p is the density matrix of the cavity state and a and the usual annihilation and 
creation operators of the photons, respectively. 

For the simplicity and without loss of the generality, let us assume that the cavity is initially 
prepared in the following few-photon superposition state 

|^(0)) = C |0) + C 1 |1) + C 2 |2), (3) 

with the complex amplitudes: Co = |Co|e J?i , C\ = \C\\, and C<i = \C2\e lif - Then, with the matrix 
elements of Wigner operator: A (a, a*) = f d 2 ze^ a] ~ a *^~ z * <ya ~ a '^ /2n 2 , in the Fock representation 
[17] 

I i Ym / | 

(n\ A (a,a*)\m) = ^(2a) n - m e { - 2la ^L^- m \A\a\ 2 ), n, m = 0, 1, 2, ... (4) 

7T V n\ 
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here n > m, a = \ao\e %e , one can easily obtain the Wigner function of the initial state 

W(a o ,a* ,0) = -[|Co| 2 -|C 1 | 2 L?(4|«o| 2 ) + |C 2 | 2 L°(4|a o r)]e ( - 2|ao|2) 

7T 

+ — e^'^IQAIKI 2 cos(2# -<p + <f>) 

7T 

- — e ( - 2|a °l 2) |CiC 2 ||ao| cos(fl - ^)L}(4|a | 2 ) 

71 

+ -e ( - 2|ao|2) |CoC 1 ||a o |cos(^ + 0). (5) 

7T 

for the above superposition initial state. Above, L J n (x) is an associated Laguerre polynomial 
defined by [18] 

L (j) (x) = yvir {n + J)l - (6) 

{ ) h } (n-«)!(J + *)!«!■ (6) 
In what follows we discuss how such a state decay in a loss cavity by investigating the time- 
evolutions of the above initial Wigner function. 

A. Dissipative dynamics for the Wigner function 

We now consider how the above few-photons superposition state dissipates in a loss cavity 
without any thermal photon (i.e., (n) th = l/[exp(^z///csT) — 1] — > 0, for the present optical 
frequency photons and at the room temperature: hv/k B T 3> 1), which is described simply by the 
following master equation [19-20] 

— = —Kidnap + pd^d — 2dpd)), (7) 
at 

with k being the loss coefficient. Our discussion is based on the time-evolutions of the Wigner 
function, i.e., 

. <8> 

with 

dC (*' A *) = Tr[^e Aat " aA *] = KTr[(2dpd^ - tfap - ptfd)e^- &x *]. (9) 
Formally, Eq. (8) can be rewritten as 

A W{a,a*) = 2nW [lipat] (a,a*) - KW& &p] {a,a*) - nW [pfi ^ ] (a, a*), (10) 



dt 

where the symbol W^a, a*) is defined as 



W [x] (a,a*) = \ f d 2 \ C [x] (\,\*,t)e aX *- a *\ C [x] (X, A*) = Tr[x e Aat " aA *], (11) 
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withWM(a,a*) = W(a, a*), and (A, A*) = C(A,A*). Note that 



and 



4 /" rf 2 A e (QA *- Q * A) ^-C(A,A*) = a*W(a,a*), 

Tl 2 J OA 

4 / d 2 A e (« A *-^)Ac(A,A*) = -aW(a,a*), 

7H J OA* 

4 / d 2 Ae (aA *- Q * A) rC(A,A*) = -jLw(a,a*), 
it 2 J oa 

4 I d 2 \ (-A)e (QA *- Q * A) C(A,A*) = JLw(a,a*) 
Ti 2 J oa* 



we then have 



H^ [apat] (a,a*) = <i 2 A C [apa,, (A,A*) 

1 72> r * 1 * ^ 1 d d ad , 
= _ y A [ a a« + — „,_-__ + -_,C(A, A«)e- 



, *A „a\* —a* \ 



Similarly, 



W^ & \a,a*) = [ a +~][a* + IA] W M( a>a *) 
As a consequence, Eq. (10) reduces to 

whose solution reads [12] 



<i 2 ao 



2 i„_„„«.-Ktiai 



WW,t) = " 2ftt / ^e 1 "^- 1 ™ " n W(a , a*, 0), 
1 — er 2Kt J 7r 



For the cavity initial state |-0(O)) we substitute Eq. (5) into Eq. (17) and get 

W(a,a*,t) = fe(- 2 l Q ' 2 )[|C r-|C 1 | 2 (2e- 2K *-l)]L?[-f a 

71 1 



2e- 



2nt- 



2 , ,9 l?n<<=~ 2K *l 2 
+ ± e - 2 H 2 |C 2 | 2 (2e- M -l) 2 L°[ 



7r 1 1 v ' ZL 1 - 2e" 



2k* J 



7T 

Above, an integral formula [21] 

rl 2 7 
7T 



+ ^| W2 | e (-2|«| 2 -2^)| a |2 cos(2 ^ _ + 0) 
7T 

+ -|C , o C , i|e ( - 2|Q|2 - K ' ) |a|cos(^ + 0) 

7T 

+ — iGi^alef-^-^Vlcos^-^pGa^-lJe-^ + l]. (18) 



mm(m,n) 

(_22£3) n!m! 

- K)!(m - K)!(-xi) m+n - 



y -f? \T? L TT7 \ — Z Z^sTX"", ite(zi) < 0, (19) 



K = 

has been used and the unassociated Laguerre Polynomial L m (x, y): 

L m (x, y) = { -iJrH m>n (x, y), H m>n = ^^e^^l |™, (20) 

was introduced [22-23] with H nhn (x,y) being the generating function of two-variable Hermite 
polynomial. 



B. Time-dependent negativity of the Wigner function 

With the above time-evolution Wigner function, we next check how its negativity changes 
with the cavity loss. Fig. 1 numerically shows these changes with the effective time nt for the 
parameters: \C±\ = 1/3, |C 2 1 = y/2/2,9 = ip = tt, <fi = 0. Here, for convenience we define the 
Wigner function W(x,p, t) in the (x,p)-space with x = (a + a*)/2 and p = (a — a*)/ (2i). One 
can see: 

(i) Initially, the Wigner function shows obviously a negativity, i.e., at certain phase space points, 
W(x,p) < 0. This means that certain nonclassical effects (such as the anti-bunching of photons) 
can be revealed in this initial cavity state. 

(ii) With the cavity dissipation, the state of the cavity decays and the negativity of its time- 
dependent Wigner function vanishes gradually. This implies that the nonclassical properties pos- 
sessed initially would be weakened with the dissipation of the cavity. 



Kt = Kt=0.2 




-0.3 -0.3 



-2-1012 -2 ^1 1 2 

X X 

(C) (d') 



FIG. 1: Wigner functions versus phase space points, (x,p) (upper line) and (x,p = 0) (lower line), of the 
few-photon superposed state (3) for different decay times, i.e., nt = 0(a, a'), 0.2(6, b'), 0.35(c, c'),3(d, d'). 
Here, the parameters in |t/>(0)) are taken as: \C\\ = 1/3, | C2 1 = %/2/2, 9 = ip = tt, (f> = 0. 



(iii) After certain times, e.g., nt > 0.35 in Fig. 1(c), the values of the Wigner functions reveal 
the expected non-negativity, i.e., W(x,p) > 0. In this evolved state the decayed cavity should 
be classical and the corresponding Wigner functions could be explained as the usual probabilistic 
distributions. 



(iv) After the sufficiently-long dissipative time, the cavity state will decay to the expectable 
vacuum state or thermal state with the mean photon number being zero (i.e., n = 0). The Wigner 
function for such a dissipated final state should be a Gaussian distribution. Indeed, from Eq. (18), 
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we have 



W{a,a*, oo) = - e (- 2 H 2 )[|C | 2 + \C X \ 2 L\{Q) + |C 2 | 2 L°(0)] 

7T 

= l e (-2H 2 )[| Co |2 + Idl 3 + |C 2 | a ] 
7T 

= -e^H 2 ). (21) 

7T 



HI. DISSIPATIVE-DEPENDENT QUANTUM STATISTICAL PROPERTIES OF THE FEW- 
PHOTONS CAVITY INITIAL STATE 

Various nonclassical effects, e.g., squeezings on quantum fluctuations and sub-Poisson photon 
statistics, in quantum optical states have attracted considerable and continuing interests [24-26]. 
Many attentions have been paid to find various non-classical optical states, while how these non- 
classical effects change with the decays of the selected non-classical states is a relatively-new 
topic. Recently, Biswas and Agarwal discussed how the Mandel Q-factor decreases with the 
decays of the photon- subtracted squeezed states[12]. Their numerical results showed that the Q- 
factor vanishes at the long dissipative times (i.e., nt — > oo) and the initial cavity state will decay to 
the vacuum. With the dissipative-dependent Wigner functions obtained in the previous section, we 
can investigate how the photonic anti-bunching effect changes with the decay of the few-photon 
superposition state |^(0) defined in Eq. (3). 

It is well-known that, if the second-order correlation function 



<? (2) (0) = (22) 



(qt 2 a 2 ) 
(ata) 2 



is less than 1, then the photonic distribution in the state is anti-bunching; otherwise, it is 
bunching. The symbol (O) represents the expectation value of the operator O in a quantum state 
p. For the present case we need to calculate the time-dependent expectation values of the operators 
a^ 2 a 2 and a)a for the decaying cavity state with time-dependent Wigner function W(a, a*, t). 
Formally, for an operator function [16] 

0(a,tf)(t) = J2 C n, m a fn (t)a m (t), (23) 

n,m 

we have 

(0(a,tf))(t) = Tr[0(a,a i )p(t)]= J d 2 aO s (a, a*)W(a, a* ,t). (24) 



On the other hand, from 

(a ] )(t) = [^ + yMA,AV)| A=A * =0 

(a)(t) = [-^-^]C(A,AV)U =A * =0 , (25) 

we can find that 

n,m 

= J rf 2 a^a, m [^ + y]"[-^-^ m e(-^^^)| A=A , =0 ^(a,a*,t) 

n,m 

d 2 aO s (a,a*)W(a,a*,t). (26) 
Comparing (24) and (26), we obtain 

O s (a,a*) = + y]"h^ - (27) 

n,m 

Specifically, if O = a) a, then 

S («,a*)lo=ata = M 2 -^ (28) 

and thus 

<**>(*) = /^(a,«-,*)0,(a, < 0|^ w 

= 4|C 2 | 2 e ( - 4Kt) + 2|C 2 | 2 e ( - 2Kt) [l - 2e ( - 2ft ' ) ] + IdlV" 2 ^, (29) 
Also, if 6 = a t2 a 2 , then 

O s («, a*)| 6=at2a2 = i - 2|«| 2 + |«| 4 , (30) 

and thus 

(a t2 a 2 )(t) = J d 2 aW(a,a*,t)O s (a,a*)\ 6=;iU 

= 2|C 2 | 2 e ( - 4K ' ) . (31) 

Above, the dissipative-dependent Wigner function shown in Eq. (18) was used. Consequently, we 
have 

(at 2 a 2 )(t) 
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(2) (0;t) 



[<ata)(t)P 



2|C 2 


2g(-4Ki) 


{4|C 2 


2 e (-4Kt) + 2 |C 2 


2g(-2K 


f )[l -2e(" 2K *)] + |Ci| 2 e(- 2Ki )} 2 



2|C 2 



[|C 1 | 2 + 2|C 2 | 2 ] 5 



/ (2) (0). (32) 
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1 2 3 1 2 3 

Kt Kt 

(a) (b) 

FIG. 2: (a): Normal-ordered correlation function g^ 2 \t) is unchanged with the decay of the few-photons 
cavity state; (b): Anti-normally-order correlation function g^ 2A \t) versus the effective decay time of the 
cavity. Here, the relevant parameters are taken as: 9 = <p = it, <j) = 0, and \d\ = VE/6, |C 2 | = y/E/3 (blue 
line), |d | = 2/9, \C 2 \ = 2/3 (red line), |Ci| = 1/3, |C 2 | = 1/3 (gray line), and |Ci| = 1/5, |C 2 | = 1/3 
(green line), respectively. 

This indicates that the normally-order correlation function g( 2 \0;t) is cavity-loss-invariant; its 
value depends only on the initial cavity state!. This is a surprise argument; imagining that the pho- 
tons in the initial cavity state is anti-bunching (i.e., g™(Q] t) < 1), then such a non-classical fea- 
ture is kept unchanged even the state approached finally to the vacuum with non-negative Wigner 
function. This argument is verified numerically by Fig. 2(a), which really shows that the value of 
<7^(0; t) is really unchanged with the decay. It is noted that, at the exact vacuum |0) the expected 
value of operator (aJa) is zero and thus the definition of g( 2 )(°) for this state is bizarre and insignif- 
icant. Therefore, our discussion always works for the dissipative process approaching to (but not 
arriving at) the exact vacuum. 

The dissipative-independence of the normally-correlation function can also be proven ana- 
lytically from the master equation (7). In fact, at any time t we have 

(tfa)(t) = Tr[a ] ap(t)} 

= (0\tfap(t)\0) + (l\tfap(t)\l) + (2|a t ap(t)|2) + ... + (n\tfap(t)\n) + ... 
= + (l|p(t)|l) + 2(2|p(t)|2) + ... + n(n\p(t)\n) + ... 

= Pll(t) + 2p 2 2^) + ••• + (*) + 

oo 

= y^np„, n (t), (33) 

n=0 
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and 



(a^ 2 a 2 )(t) = Tr[a^ a^aapit)] = Tr[a) (aaJ — l)ap(t)] 

= [(0\a^ atf ap(t)\0} + {l\a ] aa ] ap(t)\l) + ... + (n\a ] aa ] ap(t)\n) + ...] - (a f a)(t) 
= [ Pll (t) + 2 2 p 22 (t) + ... + n 2 p nn (t) + ...]- (a!a)(t) 

oo 

= ^n(n -!)/>„,„(*), (34) 



n=0 



and thus 



9( , (0;t) = m 

[(ata)(t)] 2 
(n 2 )(t)-(n)(t) 
K"X*)] a 

[pn(t) + 2 2 p 22 (t) + ... + n 2 p nn (t) + ...]- [pn(t) + 2p 22 (£) + ... + np nn (t) + ...] 
[pi(t) + 2p 22 (t) + ...+ (t) + ...}■> 

J2 {n 2 - n)p nn (t) 

= 2^5 . (35) 

n=0 

Above, p n ,n(t) is the diagonal elements of the density matrix p(t) in the Fock space. For the loss 
cavity initially prepared in the few-photon superposition state (3), one can easily see that p rhn = 0, 
for n > 2, and the other non-zero diagonal elements are determined by the following equation 
(from Eq. (7)), 

poo(t) = 2np n (t), 

pn(t) = -2K Pll (t) +4«p 22 (t), 

p 22 (t) = -4 K p 22 {t). (36) 
The solutions to these equations are 

Pn(t) = [pn(0) + 2p 22 (0)]e~ 2Ki - 2p 22 (0) e - 4K * 

P22(t)=p 22 (0)e~ iKt . (37) 

Consequently, 

2p 22 (t) 
[ Pll (t) + 2p 22 (t)f 

2p 22 {0)e- 4Kt 
[pii(0)e^* + 2p 22 (0)e- 2K *] 2 

9 {2) (0)- (38) 



2 P22(0) _ (2) 



[Pii(o) + 2p 22 (0)] 
11 



Suppose that any non-classical effect should vanish due to the dissipation, the dissipative- 
independence of the normally-correlation function implies that such a parameter should not be 
a good physical quantity to describe the cavity loss. Alternatively, the anti-normal ordered corre- 
lation function, defined as 

<2 A) rm _ (S^> _ (a t2 a 2 )+4(atq) + 2 
9 {) (aat)2 [(ata) + l]2 ' 

could be utilized to describe the dissipative process of the few-photon cavity. Indeed, with 
Eqs. (29) and (31) we have 

{2A) _ (aVa 2 )(t)+4{tfa)(t) + 2 
9 lU ' tj " [(ata)(t) + ip 



2\C 2 


2 e (-4«t) _ 


h4{4|C 2 


V 


-4«t) + 2 \C 2 


2 e (-2«t)[ 1 _ 2e (-2*t)] + 




2 e (-2^)} + 2 


{4|C 2 


2 e (- 4 "*) + 2|C 2 


2 e (-2«t)[i _ 2 e (" 2K *)] + 




2g(-2«t) _|_ ]_p 



4|Ci 


2 e (-2«t) + g 


c 2 


2 e(' 


-2«t) + 2 \c 2 


2 e (-4«t) + 2 


[\Ci 


2 e (-2«t) +2 |C 2 


2g(-2fct) _|_ ]_j2 



(40) 

which is not an invariant during the cavity dissipation. One can see also from Fig. 2(b) that, the 
value of the anti-normal correlation function changes with the cavity loss. After a sufficiently-long 
decay time the value of g^ 2A \0; t) should limit to 2, whatever its initial value is less than 2 or not. 
Certainly, such a dissipative-dependent behavior of the g( 2A \0; t)-parameter can also be exactly 
verified by using the analytic solutions, i.e., Eq. (37). In fact, we can see that 

(aA) (at 2 a 2 )(t)+4(atq)(t) + 2 

9 [<ata)(t ) + l] 2 



4pn(0)e- 2Kt + 8p 22 (0)e- 2Kf + 2p 22 (0) e - 4Kf + 2 
[pn(0)e- 2Kt + 2p 22 (0)e- 2K * + l] 2 

It is consistent with the Eq. (41), as if t — > oo, Eq.(42) can be shown 



(41) 



g {2A \0;t^ oo) =2. (42) 



IV. POSSIBLE EXPERIMENTAL VERIFICATION: THE PREPARATION OF FEW-PHOTON 
SUPERPOSED STATES AND MEASUREMENT OF ITS WIGNER FUNCTION 

We now discuss how to test the above arguments with a typical cavity QED system, i.e., highly 
excited Rydberg atoms in a high-Q microwave cavity [28]. An ideal setup is schematized in Fig. 4, 
wherein an atom is emitted from the source O and then flies across sequentially a quantized cavity, 
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FIG. 3: An experimental setup for preparing the superposition states of |0), |1) and \2). Here, an atom is 
emitted from the source O, then it flies sequentially across the J-C cavity, the classical microwave field, and 
at last is detected in the detector I. 

a classical microwave field, and finally is detected in the detector /. When the atom passes through 
the quantized cavity, the usual Jaynes-Cummings model with the Hamiltonian (h = 1) 

H = uj a S z + u c d J a + g(aS + + a + S^), (43) 

works. Here, tu a , co c are the atomic transition frequency and the cavity field frequency, respec- 
tively. S z , S± are the atomic operators, such that [S+, £_] = 2S Z , [S z , S±] = ±S±. a and a) are 
the annihilation and creation operators of the cavity field, respectively. And, g is the atom-field 
coupling strength. 

Initially, the atom is in the ground state |ei) and the cavity mode in the vacuum state, i.e., the 
wave function of the atom-cavity system is \4'(0)) = |0, ei). Next, the atom is injected into the 
cavity and the state of the atom-cavity system evolves to 

\ifj(t)) 1 = cos{gt 1 )\0,e 1 ) -ism(gt 1 )\l,g 1 ), (44) 

after the passage time t\. Then, we let the atom continuously across a classical microwave field 
for evolving the atomic states as: \e\) — > cos(0i/2)|ei) + ie~ tVl sm(9i/2)\gi) and \gi) — > 
cos(0i/2)|<7i) + ie lipi sin(#i/2)|ei). Here, the values of 9\ and ipi are adjustable. Therefore, 
before arriving at the atomic detector I, the state of the atom-cavity system reads 

= [cos(^ 1 )cos(|)|0) + e^sm(^ 1 )sin(|)|l)]|e 1 ) 

+ [ie-^cos^Osin^lO) - zsm(^0 cos(^)|l)]|^). (45) 

In order to generate the desirable superposition of the states |0), |1) and |2), we must let another 
atom (as the same of the former one) pass sequentially across the cavity and the microwave field. 
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Finally, the state of the whole system including two atoms and a cavity mode can be expressed as: 

9 9 

|V , (*)>2 = |0){cos(5fti)cos(5ft 2 )cosy cosy|ei)|e 2 ) 

9 9 

+ ie~ lV2 cos(gtt) cos(gt 2 ) cos — sin y]|ei)|^ 2 )} 

9 9 

+ 1 1} {sin(^ti) cos(gt 2 ) sin — cos ^e ivPl |ei) |e 2 ) 

+ ze" ivi e^ 2 Bin^i) cos(#t 2 ) sin y sin -l\ ei )\g 2 ) 

9 9 

- i cos(gtx) sin(gt 2 ) cos — cos —\e x )\g 2 ) 

9 9 

+ e llfi2 cos(^ti) sin(#t 2 ) cos — sin — |ei) |e 2 )} 

+ |2){e^e^ 2 sin(^i) sin(^ 2 ) sin ^ sin y | ei )|e 2 ) 

- ze^ 1 sin(pti) sm(gt 2 ) sin — cos — |ei)|^ 2 }}, (46) 

As a consequence, the desirable few-photon superposed state can be generated by the state- 
selective measurements on the atoms. For example, if the atoms are detected at the state |ei) |e 2 ), 
then the cavity mode collapses into 

1 9 9 

\fp(t))2 = -/={|0)[cOs(#ti)cOs(#t 2 )cOSy COSy] 

9 9 9 9 

+ 1 1) [e ltpi sin(gti) cos(gt 2 ) sin — cos — + e tip2 cos(gti) sin(gt 2 ) cos — sin — ] 

+ |2)[e^e^ sin^) sin(^ 2 ) sin | sin |]}| ei )|e 2 ), (47) 

with 

at r @l 9 2 -.Q r 9\ 9 2l2 

N = [cos gt\ cos gt 2 cos — cos —\ + [sin gti sin gt 2 sin — sin —\ 

r ^1 ^2i9 r @1 @2-i 2 

+ [sin gtx cos gt 2 sin — cos —J + [cos gti sin gt 2 cos — sm —\ 

+ - cos(v?i — <f 2 ) sin 1gt\ sin 2gt 2 si n $1 srn $2 (48) 

8 

being the normalized coefficient. If the relevant parameters are set properly as: if\ = ir,<f 2 — 
0, gt\ = gt 2 = 9 2 /2 = n/A,9 1 /2 = 7n/A, then a typical few -photon state discussed above 

\o) + fm+^ |2> (49) 



can be obtained. 
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The method to measure the Wigner function for a given cavity state is relative standard. Here, 
we recommend the approach proposed by Lutterbach and Davidovich [28] by directly detecting the 
the negativity of Wigner function via the atomic Ramsey interfero merries. In fact, at a phase space 
point a, Wigner function for the cavity state with the density matrix p can be simply expressed by 
[29] 

W(a) = 2Tr[D(-a)pD(a)P] = 2(P). (50) 

Here, P = exp(i7rd + a) and D(a) = exp(aa + — a* a). Furthermore, the quantity (P) can be 
determined by measuring the probability P e (or P g ) of the atom is detected at its excited state \e) 
(or \g)), i.e., 

P e (0,a) = i[l + (P)cos# (51) 
Therefore, the Wigner function is determined by 

W{a) = 2[P e (0,a)-P e (ir,a)}. (52) 

Consequently, if we have 

P e (0,a) < P e (ir,a), (53) 

then the Wigner function attains a negative value. With these preparations and measurements, the 
dissipative dynamics presented above could be tested experimentally. 



V. DISCUSSIONS AND CONCLUSIONS 



With the few-photon superposed state, in this paper we have investigated the dissipative dy- 
namics of the quantized mode without any thermal photon. By numerical method, we discuss how 
the Wiginer function of the cavity state changes with the dissipation of the cavity. Our results 
show clearly that the initial negativity of the Wigner function vanishes with the cavity dissipa- 
tion. With the dissipative-dependent Wigner function, we further discuss how a typical quantum 
statistical property, the second-order correlation function g( 2 \0), changes with the dissipation of 
the cavity. It is surprised that such a quantity is an invariant during the dissipation of the cavity. 
This argument was also verified by analytical method directly solving the master equation of the 
dissipative cavity. This implies that the g^(0) should not be a good physical quantity to describe 
the dissipative dynamics of the cavity, at least for the few-photon state. 
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The discussion in the present work is limited to the photons in optical cavity, and thus the mean 
thermal photons at room temperature can be really negligible. This implies that the final state of 
the dissipative optical cavity is exactly vacuum, at which the standard definition of the second- 
order correlation function is bizarre and insignificant. The generalization to the dissipative cavity 
with non-zero thermal photons is in progress. 

Given the few-photon superposed state of the cavity is not difficult to be prepared and its rele- 
vant Wigner function can also be easily measured in the usual cavity QED system, we believe that 
our arguments are testable with the current experimental technique. 
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